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Small metallic grains which satisfy the conditions of the universal Hamiltonian are considered. It 
is shown that for such grains the effects of the interactions in the spin channel and in the Cooper 
channel on their spin magnetization are well separated, thus allowing the determination of the 
interaction parameters within this model. In particular, the existence of pairing correlations in 
small grains and the sign of the interaction in the Cooper channel can be uniquely determined. 



I. INTRODUCTION 

In general, the problem of disorder and interaction in 
electron systems is a very difficult one. However, it was 
shown [1-4] that for small diffusive metallic grains with 
large dimensional conductance g — E^h/d the problem 
simplifies considerably. Here d is the mean level spacing 
and £rh = TiD / L? is the Thouless energy which is the 
inverse time to diffuse across the grain. D is the diffu- 
sion constant and L is the grain's size. The low energy 
physics of such small grains is described to leading order 
in 1 1 g by the "universal Hamiltonian" [4] , in which only 
the diagonal matrix elements of the interaction survive: 

n/d 

H = Y,J2 e ncl <a c n>a + E C N 2 + J c ftf + J S S 2 . (1) 



The index n spans a shell of 0,/d doubly degenerate 
time reversed states of energy e n , N = J2n=i J2a c n,a c n,a 



a,a' C n,a a cr,cr' C n,a / 



is the number operator, S = | J2n=i S 

is the total spin operator, and T = J2n=i C n,-C n .+ is the 
pair annihilation operator. E c is the charging energy and 



J, 



is) 



A c ( s )<i, where A c and X s are the dimensionless in- 



teraction parameters in the Cooper channel and in the 
spin channel respectively, fl is of the order of £rh , and 
we take Q/d = 2g. Recently, a similar problem of a ballis- 
tic grain with chaotic boundary conditions was addressed 
using renormalization group approach, and it was shown 
[5,6] that for weak interactions the low energy physics is 
indeed controlled by the universal Hamiltonian. 

This relatively simple description of the low energy 
physics of diffusive metallic grains provides the opportu- 
nity to consider theoretically, and eventually experimen- 
tally, problems which in bulk systems are much harder 
to attack. One interesting problem is the question of 
whether metals such as Gold, Copper, and Silver are su- 
perconducting or not at very low temperatures [7], i.e. if 
their effective interaction in the Cooper channel is attrac- 
tive or repulsive. While all these metals are not found 
to be superconducting down to currently accessible tem- 
peratures, it may well be that their effective electron- 
electron interaction is attractive but small. Since T c 
depends exponentially on the interaction, such weak in- 
teraction will lead to unmeasurable T c . However, small 



effective attractive interaction in such metals would af- 
fect other properties, like the proximity effect [8,7] and 
persistent currents [9,10], which depend linearly on the 
interaction. Furthermore, the magnitude of the effective 
attractive interaction in these metals may be size depen- 
dent, as can be inferred from the apparent size depen- 
dence of T c in many superconducting materials [11-14]. 
In particular, Platinum, which is not known to be a su- 
perconductor in bulk form, was recently reported to be 
superconducting at very low temperatures in granular 
form [15]. 

While the determination of the effective interaction in 
bulk materials is a difficult task, it was already recog- 
nized that weak pairing correlations can be detected in 
small "superconducting" grains [16-18]. In these works it 
was shown that the existence of weak pairing correlations 
will result in measurable effects in the spin susceptibility 
[16,17] and specific heat [18] of the grains. All these works 
considered the reduced BCS Hamiltonian, in which only 
the pairing interaction exists. However, in a real system 
other interactions exist, and in order to experimentally 
determine the existence of pairing correlations one has to 
show that the measured effect is uniquely caused by the 
pairing interaction itself. 

Small disordered metallic grains with j > 1 and not 
too strong interactions [5,6] are favorable from this point 
of view, as they satisfy the validity conditions of the uni- 
versal Hamiltonian model, and therefore the constraints 
this model dictates on the interaction terms. In this pa- 
per we calculate the ensemble averaged differential spin 
susceptibility Xs at T = of such isolated grains, and 
show that the effects of the different interaction terms are 
well separated, thus allowing an unequivocal determina- 
tion of the existence of pairing correlations in such grains, 
and furthermore, a determination of the sign and magni- 
tude of the effective interaction constants as they appear 
in the universal Hamiltonian. Actually, we consider the 
determination of A s and A c only. Since the grains are iso- 
lated, the charging energy E c is not relevant, and could 
be determined by complementary tunneling experiments. 
We consider the regime of |A C |, |A S | <C 1. Note, that for 
A c < two regimes exist, the perturbative regime and the 
superconducting regime, for which |A C | > l/ln[i?Th/d] 
[17]. We first consider the former, and then the latter 
regime. 
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II. THE PERTURB ATIVE REGIME 

Using the universal Hamiltonian, wc assume that the 
spin-orbit interaction is small and neglect it [4]. This 
assumption should be verified when comparing our re- 
sults with experiments, keeping in mind the specifics of 
spin-orbit interaction in small grains (see e.g. [19,20]). 
Throughout the paper we will be interested in the ensem- 
ble averaged differential spin susceptibility at magnetic 
fields H ^> d/fiB- In this regime we can neglect level 
statistics and assume that the energy levels in the grains 
are equally spaced. Differences between grains with odd 
number and even number of electrons can be neglected in 
this regime as well, and for simplicity we consider grains 
with even number of electrons. For detailed considera- 
tions regarding the neglect of level statistics and even-odd 
effects see section III of Rcf. [17]. In particular, ensem- 
bles of the order of 10 6 grains or larger are required for 
the shift in the magnetization (see below) to be larger 
than the fluctuations due to level statistics. We also 
neglect orbital magnetization. This can be achieved in 
pancake shaped grains (see e.g. Ref. [21]), when the field 
is applied in the direction of the thin part. Practically, 
orbital magnetization can not be completely avoided, but 
its relative magnitude can be experimentally determined 
by changing the direction of the applied magnetic field. 

The spin magnetization of a grain is given by 



M = ns(n+ — n_) 



(2) 



where n+ and n_ are the number of electrons with spin 
parallel and anti-parallel to the magnetic field, respec- 
tively. We define I as the number of flipped spins, such 
that n + — n_ = 21. It can be shown that among all states 
with I flipped spins, the one that has the lowest energy 
has all I states above Ep and I states below Ep singly 
occupied by electrons with spin parallel to the magnetic 
field. The number I that is realized at a given magnetic 
field is the one minimizing the total energy of the grain: 



E{l) = E + El ri + El t -2l m H . 



(3) 



Here Eq is the energy of the noninteracting Fermi state 
(with I = 0, no singly occupied single particle states), 
EL n = l 2 d is the kinetic energy cost of flipping I pairs, 
E( nt is the energy due to the interaction, and —21{irH 
is the Zeeman energy. In order to calculate E l int we use 
Richardson's exact solution [22,23]. Although this solu- 
tion was derived for the reduced BCS Hamiltonian, it can 
be easily generalized to solve the universal Hamiltonian 
for isolated grains. The N 2 term is then not relevant, 
and the only relevant extra term in the universal Hamil- 
tonian compared to the reduced BCS Hamiltonian is the 
spin term. 

Given I flipped spins, levels g — I + l...g + I = B are 
singly occupied, and do not participate in the pairing in- 
teraction [24]. Denoting U = fi \ B, and neglecting the 
spin term, Richardson's solution is given by a set of k 



coupled nonlinear equations, the i/'th equation of which 
is given by [23]: 



1 

X r d 



E 



— E,, — E,, 



Here k is half the number of the "paired" electrons, and 
in our case k = g — I. The total energy of the system is 
given by 



£bcs = E £ j ■ + E E - ' 



(5) 



and the many-body wave function is also given in terms 
of the k energy parameters {E v } which solve the equa- 
tions (4) . Since the electrons participating in the pairing 
interaction have zero total spin, including the spin term 
and the Zeeman term does not change Richardson's equa- 
tions, energy parameters, and orbital wavefunction. The 
spin and Zeeman terms do change the energy of the sys- 
tem, for a given I by E s = X s dl(l + 1) and Ez — —21^-qH 
respectively. 

The total energy can therefore be written as: 

B k 

E V) = E c j ■ + E E » + Xsdl{ - 1 + v >~ 2l ^ H > ( 6 ) 

j v=l 

or, in accordance with Eq. (3) 

k 

E(l) = E a + l 2 d+Y^ SE V + X s dl{l + 1) - 2ln B H , (7) 

where 8E V = E v - 2e„. Therefore, E int = X s dl(l + 1) + 
E pair where 



k 



Epair = 8E V 



(8) 



is the energy due to the interaction in the Cooper chan- 
nel, and the problem reduces to finding E pa i T (l). In 
Ref. [17] this was done to second order in the interaction 
A c . Here we use Richardson's exact solution for the de- 
termination of -Epair (0- This formalism allows a rigorous 
inclusion of the spin term. It also allows the possibility 
to give a general expression for £7 pa ir(0> an d then obtain 
the result to second order in A c as an expansion of the 
exact result. 

Manipulating Eq. (4) one obtains [17] 



5E V 



Xrd 



1 + A c a v 



where 



l ^277 i X E 



E,, — E„ 



(9) 



(10) 
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For the lowest energy solution, we approximate 8E„ by 
5E° = X u d , where X v = — y — ^ , (11) 

1 + A c a v 

and a° = a„(A c = 0) is given by 

u k 

a o = V - V 

' v ^2i-2v ^ ii-v' 



(12) 



This approximation is exact to second order in A c , and its 
accuracy to higher orders in A c was studied in Rcf. [17]. 
-Epair can now be calculated to any order in A c by insert- 
ing expression (12) in Eq. (8). To second order in A c this 
gives 



1 9 ~ l 

E paiT (l) = \ c d(g - I) + -X 2 c dJ2^ 



• + 1 + V 



21 + 1/ 



■ (13) 



Inserting Eq. (13) into Eq. (7) and differentiating with 
respect to I we obtain an equation for I that minimizes 
E(l) 



2ld + X s d(2l + 1) - X c d + X 2 c dln 



which results in 



g_ 

211 



- 2/iB-ff = 



(14) 



M = 



H B [2fi B H/d - Xl \n[E Th /(2 m H)] + A c - A s 
1 + A S 



(15) 



In Eqs. (14) and (15), for the values inside the loga- 
rithm, we assume I <§C g and replace I with its nonin- 
tcracting value. The I that minimizes E(l) as obtained 
from Eq. (14) is given by the condition that the energy 
gain from the Zeeman term when flipping another elec- 
tron and creating 2 additional singly occupied states with 
spin up electrons is equal to the energy cost of flipping 
this electron, resulting from the kinetic energy, spin in- 
teraction and pairing interaction. The kinetic part alone 
produces the noninteracting result [xo m Eq. (16) below 
for the susceptibility]. The leading contribution of the 
spin part to the total energy is proportional to I 2 , like 
the kinetic energy, and this results in an effective renor- 
malization of the density of states. The second part of 
the spin term, as well as the leading part of the pairing 
interaction, contribute to the total energy terms which 
are linear in I, like the Zeeman term, and therefore result 
in a constant shift of the magnetization, and do not affect 
Xs ■ The field dependent correction to \ s comes from the 
higher orders of the pairing term, of which the second 
order gives the dominant contribution. This part gives 
a negative correction to the energy which is monotoni- 
cally decreasing with increasing I, therefore contributing 
a positive, field dependent contribution to Xs- 



Differentiating with respect to H we obtain the en- 
semble averaged spin susceptibility for d/fi B *C H <C 
-Etii/mb 



Xs = 



Xo 



i + x s 



1 



2fi B H 



(16) 



This is our central result. The interaction in the spin 
channel results in an H independent shift of the suscepti- 
bility by a factor of 1/(1 + A s ). This gives the possibility 
to determine A s , by e.g. the Sommcrfeld- Wilson ratio, 
that compares Xs to the linear specific heat coefficient. 
The interaction in the Cooper channel results in a 1/H 
correction to Xs- This correction is a finite size effect, 
as it is proportional to the level spacing. Moreover, this 
correction unequivocally signals the presence of pairing 
correlations in small metallic grains, as it does not result 
from the interaction in the spin channel or the charging 
energy, and all other interactions have 1/g smallness. In- 
terestingly, the 1/H correction does not depend on the 
sign of the interaction, and therefore exists for attractive 
as well as repulsive interaction in the Cooper channel. 
Thus, measuring Xs m small metallic grains at magnetic 
fields H 3> c?/mb determines the magnitude of A c , but 
not its sign. In order to obtain the sign of A c one has 
to look at M/H. Unlike the case in the susceptibility, 
where the first order term in the interaction is not field 
dependent, and therefore does not contribute, here, to 
leading order in A c 



M 



Xo 



l + X s 



1 + 



(A c - X s )d 



2fi B H 



(17) 



and the 1/H correction does depend on the sign of A c . 
Once X s is cither known or small, the sign of A c is easily 
determined. Note, that in principal the information given 
by Xs and by M/H is equivalent. However, their high 
magnetic field behavior is different, and therefore both 
the sign and magnitude of A c can be obtained. (Actually, 
both can be obtained from the behavior of M/H. How- 
ever, the susceptibility measurement is preferable for the 
determination of the magnitude of A c because it is inde- 
pendent of any other interaction. It is also a more precise 
measurement experimentally). The magnetic field range 
for which our treatment is valid is given above Eq.(16), 
and depends on the specific metallic grain, as well as its 
size and its dimcnsionlcss conductance. For example, for 
Copper grains of size 5 * 50 * 50nm 3 and g — 25 the level 
spacing is roughly 0.06K, the Thouless energy 1.5K, and 
therefore the magnetic field range would be between 0.1 
and 2.5 Tesla. 



III. THE SUPERCONDUCTING REGIME 

So far we considered the perturbative regime, which for 
attractive interaction corresponds to |A C | < 1/ In [E^/d] 
which is equivalent to d > A where A is the bulk gap 
in the mean field BCS approximation. In the crossover 



3 



regime, where d w A, the behavior of Xs changes con- 
siderably in the low magnetic field regime, (J-bH < d. 
However, the properties of \s at high magnetic field 
/ibH ^> A 2 jd are similar to those in the perturbative 
regime [17], and the interaction parameters can be simi- 
larly determined. The parameters of the universal Hamil- 
tonian can also be determined in the "BCS regime", 
where |A C | > 1/ In [E^/d} and the level spacing d <C A 
and can therefore be neglected. In this regime A c is easy 
to determine, e.g. by measuring the excitation gap. In 
order to determine A s in this regime we revisit the spin 
magnetization of the system. For A s = it is well known 
[25,26] that the spin magnetization of a superconduc- 
tor is zero below a value of H = A/(v / 2a*b), where a 
sharp step to the value of the spin magnetization of non- 
interacting electrons at the same H occurs. The area 
between the magnetization curves of the noninteracting 
and superconducting systems gives the condensation en- 
ergy, A 2 / (2d). We have already shown that finite A s 
changes the slope of the spin magnetization of noninter- 
acting electrons [see Eq. (16) with A c = 0]. Here we show 
that it also changes the value of H at which the step in 
the magnetization of a superconducting system occurs, 
as to keep the area between the magnetization curves to 
equal A 2 / (2d). Thus, one can determine A s in the su- 
perconducting regime by the magnetic field value of the 
magnetization step. This value of H is where the nor- 
mal and superconducting states have the same energy, 
i.e. when the equation 



l 2 d + JJ(l + !) + _- 2lfi B H = 



(18) 



has one solution. This occurs when I = A/ y/2d(d + J s ), 
or when 



H = 



V2w 



(19) 



The shift in the magnetic field value of the spin magne- 
tization step is a direct measure of A s in this regime. 



IV. SUMMARY 



We have thus shown that the determination of the in- 
teraction parameters in small metallic grains with not too 
large interactions can be done by measuring their ensem- 
ble averaged differential spin susceptibility. Such a mea- 
surement, done systematically as function of grain size, 
can shed light on the change of transition temperature 
with grain size in granular superconductors. Although 
our theory is valid for finite size grains, and can not di- 
rectly determine if a certain material is superconducting 
at low temperatures in bulk form, a systematic measure- 
ment of the interaction parameters as function of grain 
size can suggest the bulk behavior as well. 
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